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Abstract. We investigate ultracold and dilute bosonic atoms under strong transverse 
harmonic confinement by using a ID modified Gross-Pitaevskii equation (ID MGPE), 
which accounts for the energy dependence of the two-body scattering amplitude within 
an effective-range expansion. We study sound waves and solitons of the quasi-ID 
system comparing ID MGPE results with the ID GPE ones. We point out that, when 
the finite-size nature of the interaction is taken into account, the speed of sound and 
the density profiles of both dark and bright solitons show relevant quantitative changes 
with respect to what predicted by the standard ID GPE. 
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1. Introduction 

The Gross-Pitaevskii equation (GPE), which plays a relevant role in the study of 
Bose-Einstein condensates (BECs) made of ultracold and dilute alkali-metal atoms, 
is based on the assumption of a zero-range inter-atomic potential pQ. Recently, 
several experiments [2j employing the Fano-Feshbach resonance technique in cold atomic 
collisions |3j have shown that it is possible to change the magnitude and the sign of the 
scattering length a s by using an external magnetic field. Thus, by using Fano-Feshbach 
resonances it is now possible to explore, at fixed density n, regimes where the GPE and 
its assumptions lose their validity. 

In this work, going beyond the Fermi pseudopotential approximation (contact 
interaction) of the standard GPE, we focus on sound waves and solitons in a BEC of 
interacting bosons at zero temperature under a strong transverse harmonic confinement. 
We take into account the dependence on the energy of the two-body scattering amplitude 
employing the effective-range expansion illustrated by Fu and et al. in [4] by inserting 
therein the correction proposed by Collin and co-workers in [5j. These two ingredients 
allow us to write a modified version of the Gross-Pitaevskii equation (MGPE, as named 
in ([!])) which incorporates the finite-range nature of the inter-atomic interaction. We 
reduce the dimensionality of the 3D MGPE by integrating out the degrees of freedom 
in the radial plane and we obtain a ID MGPE which takes into account both the 
scattering length and the effective range of the inter-atomic potential. We model the 
boson-boson interaction by means of three potentials: hard-sphere potential, Van-der- 
Waals potential, and square-well potential. We set the s-wave scattering length to a 
given value and calculate, for this a s , the effective range of each above model potential. 
In this way, we find relevant quantitative changes of the atomic cloud properties, i.e. 
the speed of sound and the width of the dark and bright solitons, with respect to the 
results provided by the familiar one-dimensional Gross-Pitaevskii equation. 

2. The modified Gross-Pitaevskii equation 

We consider N interacting bosons of mass m confined by an external trapping potential 
V tr ap{r) at zero temperature. The Hamiltonian is then given by 

N N 

ff = £ft(n) +j £!>(?;-f,), (l) 

i= 1 i =1 j^i 

where 

Kd) = -T^V, 2 + (2) 

2m 

with V(r t — fj ) describing the interaction between two bosons at positions r) and f r 
The ground-state properties of a weakly interacting bosonic gas can be very efficiently 
described by using the standard Gross-Pitaevskii equation (GPE) p. As well known, 
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one can derive the GPE minimizing the GP energy functional Eqp- 
inter-atomic potential by the Fermi pseudopotential 

V F {fi - ff) = g8(fi - ff) , 
where the coupling strength g is 
Airh 2 a q 


9 = 


m 


Describing the 

( 3 ) 

( 4 ) 


with a s the interparticle s-wave scattering length, the energy functional Egp reads: 

E GP [(f), (f)*] — N I d 3 f <f{f)*h{f)<f{f) + ^N(N — 1) J d 3 f](j){r) | 4 , (5) 

where <fi(f) is the single-particle wave function (all the N bosons are in the same single¬ 
particle state). By exploiting the variational approach, where the functional Eqp is 
required to have a minimum with respect to <f{f) obeying the normalization condition: 


d 3 f o(r) = 1 , 


( 6 ) 


by using that for very large N one can write that (N — 1) ~ N and by employing the 
Lagrange multipliers method, one arrives to the standard GPE 

- + V tr ap(f) +gN\(j){r)\ 2 <j>(f) = /i <f>(r) , (7) 

where /i is the chemical potential. 

At this point some considerations about the inter-ato mi c potential (J3J) are in order. 
Such a potential ignores completely the dependence on the energy of the scattering 
amplitude. This approximation, however, is valid provided na 3 is sufficiently small. On 
the other hand, for stronger confinements and larger values of na 3 s) a better treatment 
of atomic interactions that preserves much of the structure of the GP theory is possible. 
This goal can be pursued by introducing an effective interaction potential V e ff which 
gives the energy dependence of the scattering amplitude through an effective-range 
expansion which will also depend on the effective range r e of the inter-atomic potential 
m- Specifically, in the following, we use the effective interaction potential 

Veffin - fj) = V F {fi - Tj) + Vraodifi - ff) , ( 8 ) 


where 


and 


Vmodifi - ff = y [5 (d - rf)V%._^. + - fj)] 


92 = 


AttH 2 


m 




( 9 ) 


( 10 ) 


3 2 

In this case, from Eq. (l8]h it can be deduced that the energy functional has an extra 
term E mo d , due to V mo d, having the following form: 


J mod\ 


= — [ d 3 R 


~yJ - r 2 )0(fi)0(f 2 ) = 

f d 3 ff*(R + r -)<j>*{R - r -)V mod (f)f(R + r -)<t)(R - T -) , (11) 
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where we have made use of (N—l) ~ N and the second row is a re-writing of the first one 
in the two body center-of-mass frame (f = fj — fj, R = (f + fj)/ 2). The simplification 
of E mo d achieved by doing calculations in the above frame and minimization of the 
(inclusive- E mo f) modified Gross-Pitaevskii (MGP) energy functional 


Emgp[(1>*i 0 ] 




^V 2 

4 


0 ( 12 ) 


with respect to 0* with the constraint (J6]) provide the following modified Gross-Pitaevskii 
equation (MGPE) 


- ^V 2 + V trap {f) + g N\f(f)\ 2 + y£/ 2 V 2 (|0(r)| 2 ) 0(f) = //0(f) . (13) 

Notice that a similar nonlinear Schrodinger equation has been derived and studied by 
Garcia-Ripoll, Konotop, Malomed, and Perez-Garcia [6]. Their investigation starts 
from the Hartee equation for bosons, which is a nonlocal integral Schrodinger equation 
(nonlocal GPE) [7], and it is based on a gradient expansion of the nonlocal GPE (Eld- 


3. The one-dimensional MGPE 


We assume that the external confinement potential V trap (f) is obtained by superimposing 
to a very strong isotropic harmonic confinement in the x — y (radial) plane a generic 
shallow potential along the z (axial) direction, so that 

Vtrapif) = ^muj 2 ± (x 2 + y 2 ) + U(z) , (14) 

where cn_i_ is the trapping harmonic frequency. The spatial degree of freedom in the radial 
plane is thus frozen and the system can be considered, in practice, one-dimensional (ID) 
in the axial direction. As suggested by the form (fT4j) of the external trapping potential, 
we shall use the following Gaussian ansatz for the single-particle wave function 0(f): 


0(f) 


fO) 

y/n a± 


e 


x 2 +y 

2a' 2 . 


2 


(15) 


where a± = y/KJ(muj±) is the transverse characteristic length of the ground state 
of the harmonic potential and f dz\ip{z)\~ = 1. This ansatz will be valid when 
g|<^| 2 /27ra^ <C 2haq_ |8j. Inserting Eqs. (1T41) and ffT5l) into Eq. (fT2l) and then minimizing 
with respect to ip* leads to the ID version of the modified Gross-Pitaveskii equation 


h 2 d 2 
2m dz 2 


U{z) 


7 + 


where 


7 = 


2-TTa i 


9 ~ 


92 


72 = 


92 


2ira 2 , 


<p(z) = fup(z) , 


jl — fi — fruj_ j_ 


(16) 


(17) 


The effective-range effects heralded by Eq. (TT6l) become clear when the ratio of the 
absolute value of the effective range \r e \ to the inter-atomic distance (referred to the 3D 
system) is of the same order of magnitude of the ratio of this distance to the absolute 
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value of the s-wave scattering length |a s |. In this situation the dependence on the energy 
of the two-body scattering amplitude (see, for example, mini) cannot be neglected and 
the usual ID GPE is not able to describe adequately anymore the physics of our system. 
Thus, to study the effects of the finite-size nature of the boson-boson interaction on the 
atomic cloud properties, in our forthcoming ID MGPE-based studies, |r e | and |a s | will 
be chosen in such a way to meet the condition mentioned above. Moreover, note that 
results from Eq. (1161) are reliable as long as fV|a s |/a_i_ <C 1. 


4. Interaction potentials 


In this section we present three toy models for the two-body interaction potential 
between atoms. Then, we shall use these three potentials in the analysis of the sound 
velocity and solitonic waves within the system under investigation. 


• Hard-sphere potential. This model for the description of the boson-boson interaction 
is defined as follows 


V(r) = oo r < a s , while V(r) = 0 r > a s . (18) 

For this potential, 

2 . . 

= ~a s , (19) 

and one thus reduces to the standard GPE since 72 = 0, as it can be seen from see 
the first and second formula of Eq. (TT71) with $2 given by Eq. (ITOl) . 

• Square-well potential. In this case, the two-body collisions are described by a 
potential well characterized by a finite depth: 

V(r) = —V 0 r < r 0 , while V(r) = 0 r > r 0 (20) 


with Vo positive. It is possible to show that in the limit of sufficiently small incident 
wave vector (q —> 0 ), the s-wave scattering length a s is given by 

tan (x( 0 )r 0 )' 


a s = r 0 


x( 0 ) r o 


and the effective range r e by 


r e = r 0 
where x(0 ) 2 = mV 0 /K 2 . 


3 a 2 , x(0) 2 “*n> 


( 21 ) 


( 22 ) 


• Van-der-Waals potential. When the interaction is Van der Waals-like, the 
interaction potential may be approximated by a potential well for r < r 0 (this 
latter being called empty-core radius), while by a function of the form — C§/r 6 
otherwise, that is 


V[r) = 00 r < r 0 , while V(r) = —C 6 /r 6 r > r 0 , (23) 


where Cq is a parameter which quantifies the interaction strength. Note that 
the potential above is reminiscent of the Ashcroft pseudopotential used to treat 
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conduction electrons in alkali metals. For the potential (123D the s-wave scattering 
length a s and the effective range r e have the following expressions |9]: 


Cia - 


r, = 


p2 (!) 


7T 

27T 


(1 - tan <f>) l vd , 
1 + tan 2 <F , 


(24) 

3r 2 (|) (l-tan$) 2 ^’ (25) 

respectively. In the above formulas l vd is a Cg-dependent characteristic length and 
<F a function depending on the ratio l ld/ r o ; 


l nrl. 


mC(, 

~W 


1/4 


/ 2 

$ — -M. 
2 r 2 


37T 


(26) 


The forthcoming analysis will be focused on the sound velocity and solitonic density 
profiles for each of the three boson-boson interaction potential models above presented. 
We keep fixed the scattering length a s and calculate the effective interaction range r e 
by using the formulas above provided, that is, Eq. (fTOl) for the hard spheres potential 
(TT8|) . Eqs. (12T|) and ([22]) for a given Vo in the case of the square-well potential (1201) . and 
Eqs. (1241) and (!25l) for a given Cg in the case of the Van-der-Waals potential ()23lh 


5. Sound velocity 

We want to gain physical insight both in the spatial and temporal evolution of our 
system. The theoretical tool which permits us to do this is the time-dependent version 
of the modified one-dimensional Gross-Pitaevskii equation (fTfilh We suppose that 
U(z ) = 0, and scale lengths, times, and energies in units of a±, l/cv±, and hu j_, 
respectively. We use thus the following adimensional time-dependent ID MGPE: 


. d . 


1 d 2 


d 2 


2 d, 2+7M + 2 72 ^^ 


v(z,t) ■ 


(27) 


where, for simplicity of notation, we have denoted the dimensionless quantities by the 
same symbols used for those with dimensions. We are interested, in particular, in the 
consequences of a perturbation, with respect to the equilibrium, created at a given 
spatial point of the system at a given time. We start writing <p(z,t) as: 

<p(z,t) = s/n(z,t)e iS(2 ’ t} , (28) 

with n(z,t ) describing the density profile and S(z,t ) related to the velocity field v(z,t ) 
via the relation 


v(z,t) = . 


(29) 


By inserting the two equations above in the time-dependent ID MGPE (1271) . one obtains 
the hydrodynamic equations (HEs) 

2 " 


dv 

dt 

dn 

dt 


d 

dz 

d 


-v + yn + 


72 _ 4n 


d 2 1 / dn\' 

dz 2 8n 


= 0 


+ — (nv) = 0 
dz 


( 30 ) 
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At this point, let us suppose to perturb the system with respect to the equilibrium 
configuration characterized by n(z,t) = n 0 and v(z,t ) = v 0 = 0 : 

n(z, t) — no + Sn(z, t ) 


v(z, t) — v o + Sv(z, t ) . 


(31) 


We use these formulas in the hydrodynamic equations (1301) and assume to be in the 
stationary regime, Vo = 0. Under the hypothesis that the perturbation is sufficiently 
weak so as to retain only the hn-first-order terms in the HEs, we get 



(32) 


If the perturbation is a plane wave, that is 8n(z,t ) = Ae l ^ kzZ ~ ut ^ + A* e~ l ^ kzZ ~ ut \ the 
relation of dispersion which characterizes the oscillations associated to the wave induced 
by the perturbation is 



(33) 


which depends on the equilibrium density no and contains information about two-body 
collisions via 7 and 72 , see the first two formulas of Eq. ffTfjl . and Eqs. (j3J) and fflQjl . 


The perturbation will stable with respect to time for real oj that is always guaranteed 


when a s = 2/3 r e . If this is the case, the dispersion relation (|33]l is the usual Bogoliubov 
dispersion, that is 



(34) 


which, in the limit of sufficiently small wave vector (k —> 0 ) gives back the usual 
dispersion relation of the sound wave, that is 


(35) 


u = c s k 


with the velocity c s = yfn 0 r y of sound propagating in the system related to the interaction 
parameters, equilibrium density, and harmonic trap characteristics. To see more clearly 
such a dependence we use the standard units of measure so that one has 



(36) 


where we have take into account the definitions of 7 , g and g 2 . 

As above commented, we study the sound velocity c s as a function of the equilibrium 
density no, Eq. (l36j) . and analyze such a quantity for each of three interaction potentials 
previously described. 

Fig. 1 shows the sound velocity c s as a function of the axial equilibrium density no 
on varying the shape of the inter-atomic interaction potential, see Sec. 4. We have fixed 
the s-wave scattering length a s and calculated [given r 0 and Vo for the potential ( 120 |) 
and Cq and r 0 for the potential (1231) ] the value of r e for each inter-atomic potential by 
using Eq. (1T9)) for the hard-sphere potential, Eqs. (|2l]) - f)22|) for the square-well potential, 
and Eqs.(l24lb(l25l) for the Van-der-Waals potential. 
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Figure 1. Sound velocity c s vs axial equilibrium density no for a s = 0.1. Solid line: 
hard-sphere potential (flSl) [this curve is the same provided by the standard ID GPE], 
Dot-dashed line: square-well potential (121)1) [ro = 0.8, Vo = 31.05]. Dashed line: Van- 
der-Waals potential (l23ll [Co = 0.07, r 0 = 0.278]. Lengths in units of aj_, times in 
units of 1 /w_l, c s in units of aj_a;_L, no in units of l/aj_, Co in units of Hu>±a e j_. 


For any chosen set of parameters of the inter-atomic potential nnder investigation 
the final result, will only depend on the obtained value of a s and r e . Clearly, except the 
case of the hard-core potential, fixing a s several parameters of the inter-atomic potential 
under investigation will give the same r e and the same sound velocity c s . 

We observe that the behavior of the sound velocity, when the type of boson-boson 
interaction changes, is qualitatively the same. However, at a given no, by increasing 
72 > 0 one gets a larger sound velocity c s . 

The solid line of Fig.l represents the sound velocity as a function of the axial 
equilibrium density when the interaction between the bosonic atoms is described by the 
hard-sphere potential (TT9|) . Since r e = 2/3 a s - Eq. (Tl9|l - y 2 = 0 (see Eq. (HO and the 
third formula of Eq. (fT71P so that one reduces to the same behavior predicted by the 
ID GPE with a Dirac-delta interaction characterized by the assigned a s , see Eq. (j27J). 

For instance, Fig. 1 compares sound velocity versus density in the three potentials of 
interest. We can thus conclude that the finite-size nature of the inter-atomic interaction 
has the effect to produce quantitative changes in the behavior of the sound velocity c s 
with respect to that predicted by the familiar ID GPE. 
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6. Solitons 

We start by considering the time-dependent ID MGPE (l27lh When 72 = 0 we reduce 
to the standard time-dependent one-dimensional Gross-Pitaevskii equation. It is well 
known that this equation admits the possibility of studying topological configurations 
of the Bose-Einstein condensate like solitonic solutions (solitary waves preserving their 
form and propagating with a constant velocity v ) with positive (repulsive inter-atomic 
interaction) or negative (attractive inter-atomic interaction) s-wave scattering length a s 

m 

<p(z, t) = f(z - . (37) 

The solutions corresponding to a s > 0 are the dark solitons. The axial density |/| 2 
of these solitons assumes the same finite value when x —> ±00 (with x = z — vt the 
comoving coordinate of the soliton) and is characterized by an hole-structure with a 
minimum at x — 0. The difference between the phases of the wave function at ±00 is 
hnite. For a s < 0 one has the bright solitons that set up when the negative inter-atomic 
energy of the BEC balances the positive kinetic energy so that the BEC is self-trapped 
in the axial direction. In this case |/| 2 goes to zero when x —>■ ±00 and exhibits a 
pulse-structure with a maximum at x — 0. The difference between the phases of the 
wave function at ±00 is zero. 

We focus on solitary waves when the the effective-range correction is taken into 
account, that is with q 2 Hnite. Proceeding thus from the ID MGPE, we look for its 
solutions of the form (j37l) which inserted in Eq. ()27l) provide the following differential 
equation: 

-\f" + 1 f+\^(f 2 )"f = nf, (38) 

< 9 2 

where " = ——. We observe (see the discussion in the sequel) that ID MGPE 
ox 2 

admits dark (bright) solitonic solutions when the nonlinearity 7 is positive (negative). 
Therefore, due to the form of 7 - first formula of Eq. (TT7|) - it is possible to have a given 
type of soliton irrespective of the sign of a s . 

6.1. Dark Solitons 

We study the black solitons that are dark solitons characterized by a vanishing axial 
density at x = 0 and zero velocity v with respect to the condensate. It is possible to 
achieve a relation which implicitly defines the solution / of the differential equation (1381) 
that reads 

y/l - 27 2 f(z) 2 arctanh{^-) = Jff^z (39) 

Jog 

with foo being the absolute value got by / at ±00 and 7 > 0. Since 0 < \f(z)\ 2 < 1, 
the dark solitons solution exists when — 00 < 72 < 1 / 2 . 
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The density profile f(z) 2 can be thus studied as a function of the axial coordinate 
z by solving numerically Eq. (139|) when one knows the features of the boson-boson 
interaction, i.e. both 7 and 72 . To set these two quantities, we have followed the same 
procedure followed to obtain Fig.l (see Sec. 5). We have thus plotted f{z) 2 versus z, 

Fig-2. 

We observe that when one takes into account the finite-size nature of the inter¬ 
atomic interaction, the width of the solitary wave under investigation is qualitatively 
the same of that one would found by using the familiar one-dimensional Gross-Pitaevskii 
equation (solid line, see the discussion in Sec. 5) but its magnitude meaningfully changes 
with respect to the latter case. 



Figure 2. Axial density profile f(z) 2 of the black soliton vs axial coordinate 2 for 
a s = 0.1. Solid line: hard-sphere potential (IT81) [this curve is the same provided by the 
standard ID GPE]. Dot-dashed line: square-well potential (EU1) [r 0 = 0.8, Vq = 31.05]. 
Dashed line: Van-der-Waals potential E51) [Cg = 0.07, r 0 = 0.278]. Lengths in units 
of aj_, energies in units of Cq in units of Swj_a[]_, f(z) 2 in arbitrary units. 


Actually, the width A 2 at half-minimum of the dark soliton can be easily calculated 
from Eq. (ITJl) setting = 1, f(z) = 1/2, and 2 : = Az/2. I 11 this way we immediately 
find 


Az 


2 

arctanh{}f) 



(40) 


Taking into account the definitions of 7 and 72 , Eq. (TT71) with Eqs. (S} and (TT0|) . this 
formula gives the width A^ of dark solitons as a function of the scattering length a s , 
effective range r e , and transverse width a_ j_ of the harmonic confinement. 
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Figure 3. Axial density profile f{z) 2 (at t = 0) of the bright soliton vs axial coordinate 
z for a s = —0.1. Solid line: hard-sphere potential (fTSl) [this curve is the same provided 
by the standard ID GPE]. Dot-dashed line: square-well potential (121)1) [ro = 0.5, 
Vo = 82.1011]. Dashed line: Van-der-Waals potential (l23l) [C§ = 0.07, ro = 0.2492], 
Lengths in units of aj_, energies in units of Hu>±, Cq in units of Huj±a 6 L , f(z) 2 in 
arbitrary units. 


6.2. Bright Solitons 


We start from Eq. (J38]). When 7 < 0, the constant of motion for this equation 


a- = h/T+ nr- - V - k/7 


(41) 


By requiring that / and its first derivative tend to zero at ± 00 , we get K = 0. By 
imposing that / is maximum for x = 0, we obtain /1 = — y|7|/(0) 2 , and by defining 
/ = 0(x ) 1//2 we get, from Eq. (HI]) . 



8 (K - /if + ^70 2 ) 



Then, by integrating the above expression with + and by using K 
— 1 / 27 |/( 0 )| 2 , one has that 


(42) 


0 and n = 



dy* 


272 y 


y 2 {f(0) 2 -y) ' 


(43) 


The integral at the right-hand side of Eq. (H3j) can be numerically solved by allowing 
for a study of the density profile f(z ) 2 of the soliton as a function of the axial coordinate 
z setting both 7 and 72 . Therefore for the bright solitons as well, we have studied 
the density profile f{z ) 2 as a function of the axial coordinate varying the boson-boson 
interaction potential by following the same path as for the black solitons. These results 
are enclosed in Fig. 3. From the plots therein, it can be observed - as for the sound 
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Figure 4. Width Az of dark solitons (solid line) and bright solitons (dashed line) as 
a function of the coupling 72 . We set 7=1 for dark solitons and 7 = —1 for bright 
solitons. A z in units of a±, 7 in units of hui±a ±, 72 in units of huj±a^ L . 


velocity and the dark solitons - that the width is quantitatively affected by the nature 
of the inter-atomic interaction potential. The width A z at half-maximum of the bright 
soliton can be calculated from Eq. (1431) setting /(0) = 1, f(z ) = 1/2, and z = Az/2. In 
this way we immediately find 


Az 


1 f 1 / 1 - 2722/ 

VRA/4 V 2/ 2 (! — 2/) ' 


(44) 


This formula is more complex than Eq. 11401) . but Fig. [4] shows that Eq. (|40l) has the 
same behavior of Eq. (1441) once the signs of 7 are taken into account. 


7. Conclusions 

We have considered a system of interacting atomic bosons confined in a strong harmonic 
confinement in the radial plane plus a weak potential along the axial direction at zero 
temperature. We have carried out our analysis going beyond the Fermi pseudopotential 
approximation and described the gas evolution by employing a modified one-dimensional 
Gross-Pitaevskii equation (ID MGPE) in the absence of the axial potential. By using 
the latter equation we have studied the propagation of sound waves and that of solitons 
in the system under investigation. We have used the ID MGPE to study the sound 
velocity versus the axial density and the density profiles of the solitons (black and bright) 
as function of the axial coordinate by modeling the boson-boson interaction via an 
hard-sphere potential, a square-well potential, and a Van-der-Waals potential. We have 











Effective-range signatures in quasi-ID matter waves: sound velocity and solitons 13 


performed our investigations by fixing the s-wave scattering length a s and calculating the 
effective-range r e corresponding, for this a s , to each inter-atomic potential. This analysis 
has allowed us to conclude that the effective-range signatures reflect in important 
quantitative changes (with respect to the results of the familiar ID GPE) of the speed 
of sound and solitary waves density profile. 

GM and LS acknowledge financial support from MIUR (PRIN Grant no. 
2010LLKJBX). 
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